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Abstract. We fix the supports A = (Ai, . . . , A k ) of a list of tropical poly- 
nomials and define the tropical resultant 7~TZ(A) to be the set of choices of 
coefficients such that the tropical polynomials have a common solution. We 
prove that TTZ(A) is the tropicalization of the algebraic variety of solvable 
systems and that its dimension can be computed in polynomial time. The 
tropical resultant inherits a fan structure from the secondary fan of the Cay- 
ley configuration of A and we present algorithms for the traversal of TTZ(A) in 
this structure. We also present a new algorithm for recovering a Newton poly- 
tope from the support of its tropical hypersurface. We use this to compute the 
Newton polytopc of the sparse resultant polynomial in the case when T1Z(A) 
is of codimension 1. Finally we consider the more general setting of specialized 
tropical resultants and report on experiments with our implementations. 

1. Introduction 

Wc study generalizations of the problem of computing the Newton polytope of 
the sparse resultant combinatorially, without first computing the resultant poly- 
nomial. The input is a tuple A = (A±, A2, ■ . . , Ak) of integer point configurations 
in Z n . The sparse resultant 1Z(A) of A, or the variety of solvable systems, is the 
closure in (C*) Al x (C*) A2 x • • • x (C*) Ak of the collection of tuples of polynomials 
(/l) h, ■ ■ ■ i fk) such that fx = fa = ■ ■ ■ = fk = has a solution in (C*) n and each 
fi has support A t . This variety is irreducible and defined over Q |Stu94| . If TZ(A) 
is a hypersurface then it is defined by a polynomial, unique up to scalar multiple, 
called the (sparse) resultant polynomial of A. Its Newton polytope is called the 
resultant polytope of A. 

In Rcf . [Stu94j , Sturmfels gave a combinatorial description of the resultant poly- 
tope, giving rise to a combinatorial algorithm for computing its vertices from the 
vertices of the secondary polytope of the Cayley configuration Cay (.4). A drawback 
of this construction is that the secondary polytope typically has far more vertices 
than the resultant polytope. There have been attempts to compute the resultant 
polytopes without enumerating all vertices of the secondary polytope [EFK10] . A 
main contribution of our paper is an algorithm (Section l2.5|) for traversing the trop- 
icalization of 1Z(A) as a subfan of the secondary fan of Cay (A). This approach also 
allows us to compute tropicalizations of resultant varieties of arbitrary codimension. 

The tropical resultant TTZ(A) consists of tuples of tropical polynomials having 
a common solution. We show in Theorem 12.41 that T1Z(A) coincides with the 
tropicalization of JZ(A). The tropical resultant is combinatorial in nature, and wc 
present in Theorem 12.91 a simple description of it as a union of polyhedral cones, 
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each of which is the sum of a positive orthant and a linear space. The tropicalization 
of a variety is a polyhedral fan of the same dimension as the original variety. We 
derive a new formula for the codimcnsion of the (tropical) resultant in Thcorcm l2.23l 
and show that it can be computed in polynomial time using the cardinality matroid 
intersection algorithm. 

Specialized resultants arc obtained by fixing some coefficients of f^s and con- 
sidering the collection of other coefficients giving a polynomial system solvable in 
the algebraic torus. In other words, the specialized resultants are intersections of 
sparse resultants and subspaces parallel to coordinate subspaces. When the spe- 
cialized coefficient values are generic, the tropicalization TlZs{A) of the specialized 
resultant is the stable intersection of the tropical resultant T1Z(A) with a coordi- 
nate subspace. This is a subfan of the restriction of the secondary fan of Cay(„4) 
to the subspace and can be computed by a fan traversal. The algorithms are signif- 
icantly more complex and are described in Section [3] Moreover, using the results 
from our concurrent work on tropical stable intersections jJYllj . we describe the 
specialized tropical resultant as a union of cones, each of which is the intersection 
of a coordinate subspace and the sum of a positive orthant and a linear space. 

Computation of resultants and specialized resultants, of which the implicitization 
problem is a special case, is a classical problem in commutative algebra that remains 
an active area. In the concurrent work Ref. |EFKPllj an algorithm for computing 
Newton polytopes of specialized resultant polynomials using Sturmfels' formula 
and the beneath beyond method is presented and implemented, and the work is 
therefore highly relevant for our project. While the main focus of Ref. |EFKPli] is 
the efficiency of the computation of the Newton polytopes of specialized resultant 
polynomials, our main interest has been the geometric structure of secondary fans 
which allows traversal of tropical resultants of arbitrary codimcnsion. 

The tropical description of a polytope P is a collection of cones whose union 
is the support of the codimcnsion one skeleton of the normal fan of P, with mul- 
tiplicities carrying lengths of the edges of P. That is, the union is the tropical 
hypersurface defined by P. For example, the tropical hypersurface of a zonotope is 
the union of the dual hyperplanes (zones), and the tropical hypersurface of the sec- 
ondary polytope of a point configuration contains codimension one cones spanned 
by vectors in the Gale dual. See Section 12.31 The tropical description uniquely 
identifies the polytope up to translation, and we consider it to be an equally impor- 
tant representation of a polytope as the V- and H-dcscriptions. Furthermore, the 
conversion algorithms between these representations deserve the same attention as 
other fundamental problems in convex geometry. A contribution of this paper is an 
algorithm (Algorithm 14. ip for reconstructing polytopes from their tropical descrip- 
tions. We apply the algorithm to the tropical description of resultant polytopes in 
Theorem 12.91 to recover the combinatorics of the resultant polynomial. 

All the algorithms described in this paper have been implemented in the software 
Gfan | Jenj . Computational experiments and examples are presented in Section [5] 

2. Resultants 

Let A = (A±, A2, ■ • ■ , Ak) where each Ai = {o^i, a^2, • ■ • , o-i,mi} is a multi-subset 
of Z n , and let m = m\ + m 2 + • • • + Wfc. Throughout this paper, we will assume 
that rrii > 2 for all i. Let Q 1: Q 2 , ■ ■ ■ 7 Qk be the convex hulls of Ai,A 2 , ...,Ak 
respectively. Let (C*) Ai denote the set of polynomials of the form Y^jZi c jX aij in 
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C[x 1 ,x 2 ,.-.,x n ], where each Cj is in C* := C\{0}. Let Z C U^ =1 {C*) A ' be the 
set consisting of tuples (fi, fi, ■ • ■ , fk) such that the system of equations /i = /a = 
••• = /* nas a solution in (C*) n . 

Definition 2.1. The resultant variety, or the variety of solvable systems, is the 
closure Z of Z in IlLi( C *)' 4i and is denoted 7£(.4). 

The resultant variety is usually defined as a subvariety of rii=i or its projec- 
tivization jGKZ941 IStu94j . but we chose to work in riiLi^*)^ as tropicalizations 
are most naturally defined for subvarictics of tori. 

2.1. A simple description of the tropical resultant and its multiplicities. 

The tropical semiring T = (R, 0,0) is the set of real numbers with minimum as 
tropical addition © and usual addition as tropical multiplication 0. A tropical 
(Laurent) polynomial F in n variables x = (xi,X2, ■ ■ ■ ,x n ) is a multiset of terms 
(c, a) or c x a where c € R is the coefficient and a = (ax, a%, . . . , a n ) <G is 
the exponent. We will also write F = (J)^ a ^ eF {c x a ). The support of F is the 
multiset of a's, and the Newton polytope of F is the convex hull of its support. 

The tropical solution set T(F) of a tropical polynomial F is the locus of points 
x S R™ such the the minimum is attained at least twice in the expression 

ffi (c0i°)= min (c + a\X\ + aiXi -\ h a n x n ). 

(c,a)EF 

In other words, a point x € R n is in T{F) if and only if the minimum for (1, x) ■ (c, a) 
is attained for two terms in F, which may be repeated elements. Therefore, T(F) 
is a (not necessarily pure dimensional) subcomplex of a polyhedral complex dual 
to the (marked) regular subdivision of the support of F induced by the coefficients 
c, consisting of duals of cells with at least two marked points. See Section l2~2l for 
definitions of subdivisions and marked points. 

When F contains no repeated elements, the tropical solution set coincides with 
the non-smooth locus of the piecewise-linear function from R" to R given by x 
F( x ) = (B( c ) 6 f( c 9 xCL )i which is also called a tropical hypersurface. In particular, 
if all coefficients of F are the same and if F contains no repeated elements, then 
the tropical hypersurface is the codimension one skeleton of the inner normal fan 
of the Newton polytope of F. 

Let A ~ {A\, Ai, ■ ■ ■ , Ak) be as before, and let R" 4 * denote the set of tropical 
polynomials of the form 0™! 1 (cjj x Qaij ). 

Definition 2.2. The tropical resultant T1Z{A) of A is the subset of R m , or R Al x 
R A ' 2 x • • ■ x R Ah , consisting of tuples (Fi , F2, ■ ■ ■ , Fk) such that the tropical solution 
sets of F%, Fi, ■ ■ ■ , Fk have a nonempty common intersection in R". 

We can also consider the tropical resultant as a subset of Yli=i TP" ii_1 or 
IliU R Ai /(l, 1, . . . , 1)R, but we prefer to work with M m in this paper. 

Definition 2.3. Let k be a field and / C k[a;i, . . . , x n ] an ideal. The tropical 
variety T(I) of /, or the tropicalization of V{I), is a polyhedral fan with support 

T(I) := {uj G R ra : the initial ideal in w (J) contains no monomials}. 

For uj in the relative interior of a cone C u £ T(I) we define the multiplicity as 

mult w (T(J)) := dim k (k[Z" n C^]/(in w (I))) 
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when the right hand side is finite, in particular when to is in the relative interior of 
a Grobner cone of the same dimension as T(I). 

In this definition we refer to the "constant coefficient" initial ideal as in Ref. 
BJS+07| . except that we are picking out the terms with smallest w-degree. If 
the ideal / is homogeneous, T(I) gets a fan structure from the Grobner fan of 
/. When Cuj is the smallest Grobner cone in T(I) containing ui, the initial ideal 
in w (/) is homogeneous with respect to any weight in the linear span of C w . Hence 
after multiplying each homogeneous element of in^ (I) by a Laurent monomial they 
generate an ideal (in^/)) in the Laurent polynomial ring k[Z™ n C^\- 

Theorem 2.4. The support of the tropicalization of the resultant variety 71(A) 
coincides with the tropical resultant T1Z(A). 

A consequence is that we may identify TTZ(A) with the tropicalization of 1Z(A) 
and we define its multiplicities accordingly. 

We will use incidence varieties to give a proof of Theorem 12.41 Let the incidence 
variety be 

k 

(1) W := {(/i, fit, . . . , /*, x) : Vz : /,(*) = 0} C \{{C*) A > x (C*)", 

»=i 

and let the tropical incidence variety be the set 

k 

TW := {{F u F 2 , ...,F k ,X):Vi:X £ T(F t )} C ]JR A ' x R". 

i=i 

The tropical incidence variety is the tropical prevariety [BJS + 07j defined by the 
tropicalization of the polynomials /i, /2, • • • , /fe 5 where fi is considered as a polyno- 
mial in nii -\- n variables whose support in the n variables is Ai and whose mj terms 
have indeterminate coefficients. Even if Ai contains repeated points, the support 
of fi in m, + n variables has no repeated points. 

Lemma 2.5. The polynomials fi, fc, ■ ■ ■ , fk form a tropical basis for the incidence 
variety W, i.e. the tropical incidence variety coincides with the tropicalization of 
the incidence variety. 

Proof. Let K be the field of Puiseux series in t with complex coefficients. Let 
i<2, . . . , Fk, X) be an element in the tropical prevariety, i.e. F\, F2, . . . , F^ are 
tropical polynomials with support sets A\, A%, . . . , A^, and X G W l is a tropical 
solution for each F f . Let x = (t Xl , t X2 , . . . , t Xk ) G (K*) n . Then Fi G R ra " is 
contained in the tropical hypersurface of fi(xo) considered as a polynomial in the 
indeterminate coefficients, so by Kapranov's Theorem there is a point c G (K*) mi 
with val(c) = Fi giving fi(x ) = 0. Therefore (Fi,F 2 , . . . , F k , X) can be lifted to 
the incidence variety and lies in the tropical incidence variety. □ 

A consequence of Lemma 12.51 is that we may identify the tropical incidence 
variety with the tropicalization of W and we define its multiplicities accordingly. 
The following lemma follows immediately from the definitions. 

Lemma 2.6. The tropical resultant is the projection of the tropical incidence va- 
riety onto the first factor. 

We can now prove Theorem 12.41 



COMPUTING TROPICAL RESULTANTS 



5 



Proof of Theorem \2.4\ The resultant variety 71(A) is obtained from the incidence 
variety W by projecting onto the first factor rL^O^*) an< ^ taking the closure. 
This proves the first of the following equalities. 

T(K(A)) = T(riWj) = ir(T(W)) = tt(TW) = TTZ(A) 

The second follows from Ref. |ST08j which says that the tropicalization of the 
closure of a projection of W is the projection of the tropicalization of W. The third 
is Lemma \2. 5 1 and the last Lemma 12751 □ 

For each i = 1, 2, . . . , k let Pj be the Newton polytope of fi in W ai x R™, which 
is in turn embedded in M™ x IR ra . Each Pj is a simplex; in particular, the exponent 
of every term in /j is a vertex, and there is an edge between every pair of vertices. 
Every cone in the tropical incidence variety is a transverse intersection of normal 
cones to edges of Pi, P 2 , . . . , Pfc. 

The tropical incidence variety is 

(2) (J ((V(^)j 

(Ei,E 2 ,-,E k ) \i=l / 

where the union runs over all choices of pairs Pj of points from Ai and A/"(Pj) 
denotes the inner normal cone of the corresponding edge Pj in Pj. 

Lemma 2.7. Every maximal cone in the tropical incidence variety TW = T(W) 
has multiplicity one. 

Proof. Since every vertex of every Pj has its own coordinate, the dimension of a face 
of the Minkowski sum Pi+J^+- ■ -+Pk minimizing an w is the sum of the dimensions 
of the faces of each Pj with respect to to. The dimension of the incidence variety is 
m + n — k and therefore, for a generic u e T(W), the face of Pi + Pi + ■ ■ ■ + Pfc 
minimizing lo has dimension k and must be a zonotope. Consequently the forms 
in^/i), in w (/2), . . . , in^/fc) are binomials, each with an associated edge vector m £ 
Z m+ ". The vectors t>i,i>2, • • • ,Vk generate C^j and after multiplying each in w (/j) by 
a monomial it ends up in (in w (P)) C C[Z m+ " n C^]. Hence using the binomials to 
rewrite modulo (in^/)) we get that dimc(C[Z m+Tl n C^]/(in ilJ (/))) is bounded by 
the index of the sublattice generated by V\ , i>2, ■ ■ ■ , v n in Z m+n D . If we write the 
edge vectors as columns of a matrix, then the matrix contains a full-rank identity 
submatrix, so the sublattice has index one. □ 

Let 7r be the projection from IR m x R n , where the incidence variety lies, to the 
first factor R m . The tropical resultant is the projection of the tropical incidence 
variety, so 

T1l{A)= |J Trffl^))- 

(E u E 2 ,...,E k ) \»=1 / 

The Cayley configuration C&y(A) of a tuple A = (Ai, A2, . . . ,Ak) of point con- 
figurations in Z" is defined to be A x x {ei} U • • • U A k x {e k } in Z n x 7L k . We 
will also use Cay(,4) to denote a matrix whose columns are points in the Cayley 
configuration. 
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Lemma 2.8. Let E = (Ei, £?2, • ■ • , Ek) be a tuple of pairs from A\, A2, . . . , Ak 

respectively. Then the the following cones coincide: 

TT f P| Af(Ei)j = E> {ey' : o i3 - £ + rowspace(Cay(.4)). 

Proof. Let £7 be fixed. The left hand side consists of tuples of tropical polynomi- 
als (Fi, F2, . . . , Fk) £ Ili=i ^- Ai f° r which there is a point w £ R™ attaining the 
minimum for Fi at Ei for every i. 

On the other hand, the cone R>o{e;j : ay £ E{\ + R^Zj-ey consists of all 
F = (Fi, F2, ■ ■ ■ , Fk) such that the minimum for Fi evaluated at the point £ R n 
is attained at E{ for every i. For w £ R™ and F £ R A , 

F(x — w) = min c + a • (x — w) = (F — w ■ A)(x), 

where A denotes the matrix whose columns are points in ^4, so 

T{F) +w = T{F -wA). 

Moreover, the tropical solution set remains the same if coefficients of F are changed 
by a tropical scalar multiple. Therefore, changing the coefficients of Fx, F2, . . . , F^ 
simultaneously by an element in the row space of Cay(^4) has the effect of translating 
all the tropical solution sets together, so the set on the right hand side R>o{ejj : 
a.y ^ Ei} + rowspace(Cay(„4)) consists of all tuples (F±, F2, . . . , Fk) having a point 
w £ R™ achieving the minimum for Fi at Ei for every i. □ 

The following result gives a simple description of the tropical resultant as a union 
of cones with multiplicities. 

Theorem 2.9. The tropical resultant of A is the set 

(3) TTZ(A) = |jR> {e ii : a l} (£ E,} + rowspace(Cay(^l)) 

E 

where E = (E\, E2, ■ ■ ■ , Ek) and each Ei consists of two elements in Ai. The 
multiplicity of the cone associated to E is the index of the lattice spanned by the 
rows of Cay(-E) in rowspacc(Cay(i?)) n Z m . 

The cones in © do not have to intersect each other transversely. See Exam- 
ple 02b). 

Proof. The set theoretic statement follows immediately from Lemmas 12.61 and 12.81 
Let a be the cone corresponding to E in the incidence variety, and r = 7r(cr). Us- 
ing the multiplicity formula from tropical elimination theory |ST08j , the multiplicity 
of t in the tropical resultant is the lattice index [L T : 7r(L(j)], where L T = Rr n Z m 
and L CT = Rcr D Z m+n . The lattice L CT is defined by the following equations on 
(c, x) £ Z m+n 

c ■ (eij - e ik ) + x ■ (aij - a ik ) = for {ay, a ife } = E L 

and is spanned by the integer points in the lineality space of the tropical incidence 
variety and the standard basis vectors for ay ^ Ei. The rows of the following 
matrix span the lattice points in the lineality space of the incidence variety: 

Cay(.A) 
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Hence ir(L a ) is spanned by the rows of Cay (.A) and the e^'s for dj £ E%. □ 

In Ref. [DF S07] . the tropical discriminant is described as a sum of a tropical 
linear space and an ordinary linear space. This description carries over to the 
tropical resultant when A is essential, and in particular 1Z(A) is a hypersurface. 
Our description in Theorem 12.91 is different and also works for non-essential cases 
and non-hypersurface cases. Moreover, it is simpler, and we do not need to compute 
a nontrivial tropical linear space. 

The first summand in ((3|) plus the linear span of the last k rows of Cay(„4) is 
a tropical linear space obtained as a Cartesian product of tropical hyperplanes. 
Hence Theorem 12.91 can be rephrased as follows. Let C be the matrix consisting of 
the last k rows of Cay (.4), so the kernel of C is defined by equations of the form 
c;,i + Ci.2 + • • • + Q.m; = for i = 1, 2, . . . , k. Then the tropical resultant is the set 

(4) TTZ(A) = T(ker(C)) + rowspacc [A X \A 2 \ ■ ■ ■ \A k ] . 

Example 2.10. Consider the tuple A = (Ai, A2, A3) of the following point con- 
figurations in Z 2 : 

A, ={(0,0), (0,1), (1,0)}, 

(5) ={(0,0), (1,0), (2,1)}, 

A 3 = {(0,0), (0,1), (1,2)}. 

The Cayley configuration Cay(Ai, A2, A3) consist of columns of the following 
matrix, which we also denote Cay (A): 





/0 1 1 2 





1\ 




1 1 


1 


2 


Cay(A) = 


1 1 1 










1 1 1 










^0 1 


1 




The corresponding system of polynomials consist of 






h 


= en + c u y + C13X, 
= c 2 i + c 22 x + c 2 3X 2 y, 






h 


= c 3 i + c 32 y + c 33 xy 2 . 






The point 


(0,0,0,0,1,5,0,1,5) 







is in the tropical resultant variety because the tropical hypersurfaces of the three 
tropical polynomials 

F ! = © X © Y, 
(7) =Offi(l0X)ffi(5 0X 02 0F) 

F 3 = © (1 Y) © (5 X y 02 ) 
contain the common intersection points (—1, —1) and (—2, —2). Sec FigureHJ 
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Figure 1 . A tropical hypersurfacc arrangement and its dual regu- 
lar mixed subdivision (RMS) of the Newton polytopes. The mixed 
cells are shaded. See Examples 12.101 and 12 . 18f a) . 

Consider the incidence variety defined by the ideal 

/=(/i,/2,/ 3 )CC[ C ,x ±1 ,y ±1 ]. 

The resultant variety is obtained by eliminating x and y from the system, i.e. it is 
defined by the ideal / n C[c]. In this case, the resultant variety is a hypersurfacc 
defined by the resultant polynomial 

c 12 c 23 c 31 - 2ciiC? 2 cf 3 C§ 1 C32 - cl 2 C 13 C 2 2cl 3 C 31 cl 2 + C? 1 Ci2^ 3 C3i£^ 2 - 
Cl2cf 3 C 2 lc| 3 c| 2 + CllCl2Cl3C22C2 3 c| 2 + 3c? 2 Cl3C2lcl 3 C§ 1 C 3 3 + Cl\c\ 2 C22^ C \l C ^ + 
2cf 2 Cl3ci 2 C23C3lC32C33 ~ C11C12C13C21C23C31C32C33 ~ Ci 1 Cl2C22C 2 ! 3 C3lC 3 2C 33 + 
2C12C? 3 C21C22C23C| 2 C3 3 ~ ^C\\C\ 2 C\ 3 c\ 2 C 23 c\ 2 C 33 - c\ 2 C\ 3 (? 22 C 3 xc\ 3 + 
3ci 2 Ci 3 c| 1 C 2 3C 3 lC§ 3 - CiiCi2Ci 3 C2lC22C2 3 C 3 ic| 3 - Ci 1 C 2 lC2 3 C 3 iC§ 3 - 
C12C 13 C21C| 2 C 3 2C§ 3 + CiiCi2Ci 3 C^ 2 C 3 2c| 3 + Ciicfacf^C^cfg - 
C 11 Ci 3 C2lC22C2 3 C 3 2C 3 ! 3 + cfgC^C^ - 2cncf 3 c| 1 C22c| 3 + cf 1 Cl 3 C 2 lc| 2 c| 3 . 

It is homogeneous with respect to the rows of Cay (.A). Its Newton polytope is four 
dimensional, has f-vector (15, 40, 38, 13, 1) and lies in an affinc space parallel to the 
kernel of Cay (A). 

The tropical resultant is an eight dimensional fan in R 9 with a five dimensional 
linearity space rowspace(Cay(.4)). As a subfan of the secondary fan of Cay(_4), 
it consists of 89 (out of 338) eight dimensional secondary cones, which can be 
coarsened to get the 40 normal cones dual to edges of the resultant polytope. 

2.2. Secondary fan structure and links in tropical resultants. Let A 6 

Z dxm be an integer matrix with columns oi, a 2 , . . . , a m G 7L d . We will also denote 
by A the point configuration {ai, a 2l . . . , a m }. We allow repeated points in A, as 
we consider the points to be labeled by the set {1, 2, . . . , m}, and every point gets 
a distinct label. 

Following Section 7.2A of Ref. [GKZ94 , a subdivision of A is defined as a family 
A = {Ci C A : i e 1} of subsets of A, called facets, such that 

(1) dim(conv(Ci)) = dim(conv(^4)) for each i 6 7, 

(2) conv(A) = \J ieI conv(Ci), and 

(3) for every i,j e /, the intersection of conv(Cj) and conv(Cj) is a face of 
both, and Ci D conv(Cj) = Cj n conv(Cj). 
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This notion is also called a marked subdivision by some authors, as a subdivision 
depends not only on the polyhedra conv(Ci) but also on the labeled sets Ci. The 
elements in Ci are called marked. If F is a face of conv(Cj) for some Ci G A, 
then Ci n F is called a face of the subdivision. 

For two subdivisions A and A' of A, we say that A refines A' or A' coarsens 
A if every Ci G A is contained in some Cj € A'. A subdivision is a triangula- 
tion if no proper refinement exists, and in particular, every facet contains exactly 
dim(conv(A)) + 1 elements. 

Let w : A — > R be an arbitrary real valued function on A, called a weight 
vector. We can defined a subdivision of A induced by w as follows. Consider 
the unbounded polyhedron P = conv{(a, w(a))} + R>o{e<2+i} in R d+1 , and let 
{Fi : i G 7} be its bounded facets. Then the induced subdivision is {Ci : i G 1} 
where d = {a G A : (a, w(a)) G Fi}. A subdivision A is regular or coherent if it is 
induced by some function w. The partition of R" 4 according to induced subdivisions 
is a fan, called the secondary fan of A. 

Following Section 7. ID of Ref. |GKZ94j . we can construct the secondary polytope 
of A as follows. For a triangulation T of a point configuration A, define the GKZ- 
vector 4>t G M. A as 

4>t{o) '■= V"] vol(a) 

a£T:a£a 

where the summation is over all facets a of T of which a is a vertex. 

Definition 2.11. The secondary polytope S(A) is the convex hull in R A of the 

vectors 4>t where T runs over all triangulations of A. 

Theorem 2.12. |GKZ94| § 7.1, Theorem 1.7] The vertices of 11(A) are precisely 
the vectors 4>t for which T is a regular triangulation of A. The normal fan of the 
secondary polytope T,(A) is the secondary fan of A. The normal cone of Yi(A) at 
4>t is the closure of the set of all weights w G R" 4 which induce the triangulation T. 

The link of a cone C C R m at a point v € C is 

Hnk„(C) = {u G R m | 3S > : Ve between and 5 : v + eueC}. 

The link of a fan J 7 at a point v in the support of J- is the fan 

link^) = {link^C) | v G C G F}. 

For any cone a G F, any two points in the relative interior of a give the same link 
of the fan, denoted link CT (J 7 ) . If a maximal cone C G F has an assigned multiplicity, 
we let link w (C) G link,, (J 7 ) inherit it. 

We will first show that the link of a cone in a secondary fan is a common re- 
finement of secondary fans, or cquivalcntly that a face of a secondary polytope is 
a Minkowski sum of secondary polytopes. For a sub-configuration C C A, we can 
consider the the secondary polytope of C as embedded in R" 4 by setting <j>r(a) = 
for a G A\C for every triangulation T of C. On the other hand, the secondary 
fan of C embeds in M. A with lineality space containing the coordinate directions 
corresponding to a G A\C. 

Lemma 2.13. Let A be a configuration of m points, ui G R A , and A u be the regular 
subdivision of A induced by uj. Then the face F^ of the secondary polytope of A 
supported by uj is the Minkowski sum of secondary polytopes of facets in A w . 
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Proof. Let u' £ M. A be generic and p be the vertex of the Minkowski sum picked out 
by ui 1 '. Then p = J^jP* where pi is the vertex of the secondary polytope of the zth 
facet of A w picked by (a restriction of) uj' . For each cell in A w , the vector uj + eui' 
induces the same triangulation for all e > sufficiently small. These triangulations 
together give the triangulation of A induced by uj + eu>' . For each ith facet of A u , 
Pi is the GKZ-vector of the triangulation induced by uj + suj', so their sum p is the 
GKZ vector of the triangulation of A induced by uj + eui' . Hence p is the vertex 
of in direction uj'. We can then conclude that the two polytopes are the same 
since they have the same vertex in each generic direction. 

□ 

We now define mixed subdivisions as in Rcf . [DLRSIO] . For point configurations 
j4i, A2, ■ ■ ■ , Ak in M. n , with A; = {a;,j : 1 < j < rrii}, the Minkowski sum 

k 

M = {ai.ii + a 2 ,j 2 + ■■■+ a k , jk : 1 < ji < ™i } 

2=1 

is a configuration of m\m,2 ■ ■ ■ points labeled by [mi] x [7712] x • • • x [m^]. We 
use the word cell to denote the set of labels of a face of a subdivision. 

Definition 2.14. A subset of labels is a mixed cell if it is a product of labels 
Ji x J 2 x • • • x J/j where Ji is a nonempty subset of [m<] , and it is fully mixed if in 
addition Ji contains at least two elements for every i = 1, 2, . . . , k. A subdivision 
of the Minkowski sum Ylj—i A% is mixed if every facet is labeled by a mixed cell. 

A mixed subdivision of Yli=i ^-i is also referred to as a mixed subdivision A = 
(Ai,A2, . . . , Ak). Our definition of fully mixed cell differs from that of Ref. [DFS07] 
Section 6 where it is required that conv(<Zij : j E Ji) has affine dimension at 
least one, while we only require that Ji contains at least two elements. These two 
definitions coincide if Ji contains no repeated points. 

A mixed subdivision is called regular if it is induced by a weight vector 

w : Aj — > R, where w : a,,^ w i,ji 

i—l i—1 i—1 

for some (w u w 2> . . . , w k ) E M™ 1 x K™ 2 x • • • x R mfc . In Ref. |Stu94j regular mixed 
subdivisions (RMS) were also called coherent mixed decompositions. 

Theorem 2.15. |Stu94| Theorem 5.1] For a subdivision A o/Cay(_4), the collection 
of mixed cells of the form X)»=i Q such that C\ C Aj and Uj=i C< * s a f a cet of A 
forms a mixed subdivision o/^ =1 A- TTiis is a one-to-one correspondence between 
the subdivisions o/Cay(A) and RMSs o/X)!=i ^i- Moreover the partition of weight 
vectors (wi,ui2, ■ ■ ■ ,Wk) E R mi x R'" 3 x • • • x M mfc according to the induced RMS 
coincides with the secondary fan of Cay (A). 

From our description of tropical resultants, we get the following result which 
was proven for the resultant hypersurfaces in Theorem 5.2 of Ref. [Stu94j and 
stated for the essential configurations with no repeated points in Proposition 6.8 of 
Ref. |DFS07j . See Remark 12.251 for a characterization of being essential. 

Theorem 2.16. The tropical resultant is a subfan of the secondary fan of the Cay- 
ley configuration Cay(Ai, A2, . . . ,Ak), consisting of the cones dual to subdivisions 
with fully mixed cells. 
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The multiplicities of secondary cones in the tropical resultant will be computed 
in Proposition 12.201 below. 

Proof. For a tropical polynomial F G M. A the tropical solution set T(F) is dual to 
the cells with at least two elements in the subdivision of A induced by the coefficients 
of F. More precisely, by the definition of tropical solution sets, w G T(F) if 
and only if (l,w) is an inner- normal vector for the convex hull of lifted points 
{(c, a) G M™ +1 : c x a is a term in F} supporting at least two points of A. The 
two points supported need not have distinct coordinates. 

Let (Fi,F 2 , . . . ,F k ) G R Al x R A2 x • • • x R Ak . The union of tropical solution 
sets |Jj =1 T[Fi) has a polyhedral complex structure in the common refinement of 
the completions of T(i 7 i) to W n , which is dual to the RMS of A induced by the 
coefficients of (Fi, F2, . . . ,Fk). The tuple (F±, F2, . . . , Fk) is in the tropical resultant 
if and only if the tropical solution sets have a common intersection. In other words, 
there is a fully mixed cell in the dual RMS. □ 

The tropical resultant is a subfan of the secondary fan. It is pure and connected 
in codimension one, so we can compute it by traversing. To traverse the resultant 
fan, we need to know how to find the link of a cone. 

Proposition 2.17. Let A = (A\, A2, . . . , The support of the link of a point 

lu in the tropical resultant F1Z(A) is a union of tropical resultants corresponding to 
sub- configurations of fully mixed cells in the RMS A w of A induced by oj. 

Proof. By definition, a point u is in the link if and only if uj + eu induces a RMS 
with a fully mixed cell for all sufficient small e > 0. This happens if and only if at 
least one of the fully mixed cells in A w is subdivided by u into a RMS with a fully 
mixed cell, i.e. u is in the tropical resultant of the sub-configurations of fully mixed 
cells. □ 

Example 2.18. Let A be as in Example 12.101 

(a) The link of the point (0,0,0,0,1,5,0,1,5) in the tropical resultant is a 
union of two hyperplanes whose normal vectors are: 

(0,-1,1,-1,1,0,1,-1,0) and (0,0,0,0,1,-1,0,-1,1) 

respectively. They are the resultant varieties of the sub-configurations of 
the two fully mixed cells. See Figure [TJ 

(b) The link of the point (0, 0, 0, 0, —1, —1, 0, 0, 1) consists of four rays modulo 
lineality space, three rays from the resultant of one fully mixed cell and two 
from the other, and the two resultants overlap along a ray. See Figure [2] 

The following lemma follows immediately from the definition of induced or regu- 
lar subdivisions and shows that the description of the tropical resultant as a union 
of cones in Theorem l2.9l is somewhat compatible with the secondary fan structure. 
For any tuple E — (Ei,E 2 , ■ ■ ■ ,Ek) of pairs E\ G A iy let Ce ■— R>o{ e ij : a ij £ 
Ei} + rowspacc(Cay(^4)) be the cone as in Theorem l2.9l 

Lemma 2.19. For each tuple E as above, the cone Ce is a union of secondary 
cones of C&y(A) corresponding to mixed subdivisions of Y^i=i having a mixed 
cell containing -^i- 
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FIGURE 2. The tropical solution sets at (0,0,0,0,-1,-1,0,0,1) 
and the corresponding dual RMS in Example l2.18f b). 

Let a be a secondary cone of Cay(.4) which is a maximal cone in the tropical 
resultant TTZ(A) , and let A CT be the corresponding regular mixed subdivision. Then 
all the fully mixed cells in A a are of the form Yli=i where each E is a tuple of 
pairs as above. Otherwise a would not be maximal in T1Z(A). 

Proposition 2.20. The multiplicity of the tropical resultant T1Z( A) at a secondary 
cone a of Cay (A) is the sum of multiplicities of cones Ce (given in Theorem ] 2. 9\) 
over all tuples E of pairs forming a mixed cell in the corresponding mixed subdivision 

Proof. By Lemma 12.191 for each tuple E of pairs, the cone Ce contains a if and 
only if y Ei is a mixed cell in A CT . Otherwise Ce is disjoint from the interior of 
a. The multiplicity of a is the sum of multiplicities of Cg's containing a. □ 

The edges of the resultant polytope are normal to the maximal cones in the 
tropical resultant, and Proposition ^. 20l can be used to find the lengths of the edges. 
From this description, we can derive Sturmfels' formula |Stu941 Theorem 2.1] for 
the vertices of the resultant polytope. 

2.3. Tropical description of secondary polytopes. We will give a tropical 
description of secondary polytopes of arbitrary point configurations and show how 
tropical resultants fit in. 

Proposition 2.21. Let A be a d x m integer matrix whose columns afftnely span 
an r-dimensional space. The tropical hypersurface of the secondary polytope of the 
columns of A is the set 

|J M>o{e t : i £ 1} + rowspace(A) + K{1}. 

/C{l,....m} 

|J|=r+2 

where 1 denotes the all one vector in K m . 

Proof. Let ui G K m , and let A^ be the regular subdivision of the columns of A 
induced by uj. Then ui is not in the tropical hypersurface of the secondary polytope 
if and only if A w is not a triangulation, which happens if and only if there exists a 
facet of A w containing at least r + 2 points of A. For an r + 2-subset I of {1, ... , m}, 
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Figure 3. A projective drawing of the tropical hypersurface of 
the secondary polytopc of the Cayley configuration of two 1- 
dimensional configurations in Examplc l2.22l The tropical resultant 
is shown in bold/color. A vertex labeled ij represent the vector 
dj in M 6 = R Al x M. A2 , and an edge between ij and kl represents 
the cone R>o{ey, e/y} + rowspace(Cay(.A)). Compare with dual 
pictures in Figure 2 of Ref. }Stu94j and Figure 3 of Ref. jEFKlOj . 

the cone IR>o{ej; : i £ 1} + rowspace(A) +M{1} consists of all u such that a facet of 
A w contains /. Note that for any u' <E rowspace(A) + K{1}, ui and to + u' induce 
the same subdivision of A. □ 

Comparing with Theorem 12. 9i we see that the tropical resultant is the union of 
the cones in the above tropical description of the secondary polytope of Cay (A) 
obtained by dropping an for each i. 

Example 2.22. Let A x = A 2 = {0,1,2} in Z. For A = (A 1 ,A 2 ), the tropical 
hypersurface of the secondary polytope of G&y(A) and the tropical resultant of 
A are depicted in Figure [3] The resultant polytope has f-vector (6, 11, 7, 1). The 
secondary polytopc in the case is combinatorially equivalent to the 3-dimcnsional 
associahedron and has f-vector (9, 21, 14, 1). 

2.4. Codimension of the resultant variety. In this section we discuss how to 
determine the codimension of the tropical resultant variety TTZ(A). By the Bieri- 
Grovcs Theorem BG84J this is also the codimension of 1Z(A). 

Theorem 2.23. The codimension of the tropical resultant equals 

k 



i=l 

where each Ei runs through all cardinality two subsets of A{. 

Proof. The tropical resultant variety is the collection of all lifts of all points in 
A which give a fully mixed cell in the subdivision. Therefore it is the closure 
of the collection of lifts which give a zonotope in the mixed subdivision being 
a sum of convex hull of two points from each Aj. Let P be such a zonotope 
and E — (E±, ...,Ek) the pairs of points. We wish to find the dimension of the 
(relatively open) cone Cp of lifts which will induce P. The height of the first point 
of each Ei may be chosen freely. The remaining k points of E must be lifted to the 
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same subspace of dimension dim(P), whose lift may be chosen with dim(P) degrees 
of freedom. Finally, the height of the points not in E maybe chosen generically as 
long as sufficiently large. The codimension of Cp is therefore k — dim(P). The 
theorem follows since there arc only finitely many choices for E. □ 

Lemma 2.24. Let Li denote the linear subspace affinely spanned by Ai. The 
codimension oflZA only depends on the Li and equals 

k - Max 1 , e r] i L i dim(span(v 1 , . . . , v k )). 

Proof. Since conv(£'i) C Li the quantity of the lemma is smaller than or equal to 
that of Theorem l2.23l Conversely, if we have a collection v £ Yii Li we will now show 
how we can perform a sequence of changes to v to make it only consist of vectors Vi 
which are each differences between points of Ai without lowering the dimension of 
span(w!, . . . , Wfc). Consider a vector It is a linear combination of some Uj where 
each Uj is of the form ai S — an. If all Uj belong to W := span(wi, . . . , x>i, . . . , vu) 
then so will Vi and it may be substituted by an arbitrary Uj without lowering the 
dimension. If some Uj does not belong to W then substituting uj for Vi will not 
lower the dimension. □ 

The proof also shows that instead of considering all line segments in Theorem l2.23l 
it suffices to consider only a basis for the affine span for each polytope. This is useful 
if computing the codimension using this formula. 

Remark 2.25. We can define a matroid on a set of polytopes as follows. A 
set of polytopes is independent if they contain independent line segments. It is 
straightforward to check that the base exchange axiom holds. The rank of the 
matroid is the maximal dimension of a fully mixed cell (a zonotope) spanned by two 
element subsets, one subset from each polytope. The codimension of the tropical 
resultant equals the corank of the matroid, i.e. the number of polytopes minus 
the largest dimension of such a zonotope. The (tropical) resultant variety is a 
hypersurfacc if and only if the matroid has corank one, which holds if and only if 
there is a unique circuit in the matroid. The tuple A is essential |Stu94j if and only 
if this matroid of k polytopes is uniform of rank k — 1, that is, the unique circuit 
of the matroid consists of the entire ground set. 

Using Theorem 12.231 we get a new proof of Sturmfels' formula for the codimen- 
sion of 1Z(A) . 

Theorem 2.26. [Stu94j Theorem 1.1] The codimension of the resultant variety 
1Z(A) in Y\i = i(C*y ni is the maximum of the numbers \L\ — dim(^ ?eJ Qi) where L 
runs over all subsets of {1, . . . , k}. 

By the Bieri-Groves Theorem |BG84j and Theorem 12 .41 the codimension of The- 
orem[523]cquals that of Thcorcm l2.26l In the following we explain how the equality 
of the two combinatorial quantities of Theorems 12.231 and 12 . 261 can also be seen as a 
consequence of Perfect's generalization (Theorem 12.27ft of Hall's marriage theorem 
and Rado's theorem on independent transversals. 

Let S be the ground set of a matroid with rank function p. Let U — {Si : 1 < 
i < k} be a family of subsets of S. A subset S' of S is called an independent 
partial transversal of hi if S' is independent and there exists an injection 6 : S' — ^ 
{1,2,..., k} with s G 5g( s ) for each s G S' . 
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Theorem 2.27. (Perfect's Theorem |Per69[ Theorem 2}) With the notation above, 
for every positive integer d, the family U has an independent partial transversal of 
cardinality d if and only if 

d < p(U ieI Si) + k- \I\ 

for every I C {1, 2, . . . , k}. 

In particular, the maximum cardinality of an independent partial transversal is 
equal to the minimum of the numbers on the RHS of the inequality. 

Proof of Theorem\MM Let Si = {a-b : a,b £ A,}, S = (J* =1 and U = {S. t : 1 < 
i < k}. Consider the vector matroid on S given by linear independence. Then the 
quantity MaxgdimfV,. t conv(-Ej)) is the cardinality of the maximal independent 
partial transversal of U. By Perfect's Theorem, 



Straightforward evaluation of the formulas in Theorems 12.231 and 12.261 will re- 
quire time complexity exponential in the input. Moreover, the maximal bipartite 
matching problem is a special case of this codimension problem. 

Lemma 2.28. The maximal bipartite matching problem is reducible in polynomial 
time to the problem of computing codimension of resultants. 

Proof. Let G be a bipartite graph with vertices U UV and edges E C U x V. Let 
{e u : u € U} be the standard basis for Mr . For each v £ V, let A v = {e u : (v, u) £ 
E}. Then the maximal cardinality of a bipartite matching in G is equal to the 
dimension of the resultant variety of A = ({0} U A v : v £ V). □ 

We use Theorem 12.231 to construct an efficient algorithm. 

Theorem 2.29. The codimension of the resultant can be computed in polynomial 
time in the input. 

Proof. Let A = (Ai, A2, . . . , Ak) where each Ai is a point configuration in Z". By 
Lemma \2. 241 the codimension of 1Z(A) depends only on the linear spaces L\, L2, 
. . . , Lk affinely spanned by A±, A2, . . . , Ak respectively. Choose a basis Bi for each 
linear space L^. Let B = {Bi,B2, ■ ■ ■ , B^} and 5* = U^ =1 Bi. A subset S' of S is 
called a partial transversal of B if there is an injection 9 : S' — > {1, 2, . . . , k} with 
s £ i?0( s ). The collection of partial transversals form an independent system of a 
matroid Aii on ground set S, called the transversal matroid of B. Let M.2 be the 
vector matroid on S defined by linear independence. By Theorem l2.23[ computing 
the codimension of the resultant is equivalent to computing the maximal cardinality 
of a linearly independent partial transversal, i.e. the largest subset of S which is 
independent in both M.\ and M.2- 

We can use the cardinality matroid intersection algorithm |Sch031 Section 41.2] 
to find the maximum cardinality of a set independent in two matroids with the 
same ground set. This algorithm is polynomial in the size of S and the time for 
testing independence in the matroids. Testing independence in A4\ can be reduced 
to the maximal bipartite matching problem and can be solved in polynomial time. 




Hence the two quantities from Theorems 12.231 and 12.261 are equal. 



□ 
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Testing linear independence in M.2 can be reduced to finding the rank of a matrix, 
which also takes polynomial time. □ 

The algorithm described in Theorem l2.29l is rather complex, but there is a simpler 
probabilistic or numerical algorithm. For generic vectors Vi € Li for i = 1,2, ... ,k, 
the codimension of the resultant is equal to k — rank([«i|w2 • • • \vk})- The challenge 
of turning this into a rigorous algorithm lies in making sure that the choices for 
Vi are generic. Our naive attempts at symbolic perturbations resulted in matrices 
whose ranks cannot be computed in polynomial time. 



2.5. Traversing tropical resultants. Tropical resultants are pure and connected 
in codimension 1. This allows the facets to be enumerated via the well-known 
adjacency decomposition approach. By this we mean traversing the connected 
bipartite graph encoding the facet-ridge incidences of the fan. Three operations are 
essential. We must be able to find some maximal cone in the fan, find the link at a 
ridge, and compute an adjacent maximal cone given a ray of the link at the ridge. 
In Ref. |JenlOj these subcomputations were isolated in an oracle, and we discussed 
a general algorithm for traversing a polyhedral fan (up to symmetry) represented 
only through oracle calls. In the following paragraphs we will describe how to walk 
locally in the tropical resultant. More details can be found in the next section for 
the more general setting of specialized tropical resultant. 

To find a starting cone for the traversal, we use the description of the tropical 
resultant as a union of orthants plus a linear space, as described in Theorem 12.91 
Alternatively, a generic vector in a maximal cone of a resultant fan can be found in 
polynomial time using the algorithms for the codimension described in Section |2~^1 

To find the link of a point in the tropical resultant, we use the fact that the link 
of a point w is a union of smaller tropical resultants associated to the fully mixed 
cells in the mixed subdivision of A induced by w, as shown in Proposition 12 . 1 7l 

In the tropical resultant, as a subfan of the secondary fan of Cay(.A), each cone 
can be represented by a regular subdivision of Cay(„4). The smallest secondary 
cone containing a given vector ui can be constructed from the regular subdivision 
induced by uj as explained in Section 5.2 of Ref. |DLRS10] , 

In our implementation we represent the regular subdivision A induced by ui by 
uj and the triangulation induced by a "placing" or "lexicographic" perturbation of 
lj. From this triangulation, we can easily recover the subdivision A by comparing 
the normal vectors of the facets of the triangulation lifted by w. For this to work, 
it is important to perturb ui in such a way that marked points in A remain marked 
in the refined triangulation. A full triangulation of Cay(,4) is only computed from 
scratch once at the beginning. To obtain a desired triangulation from a known 
triangulation, we find a path in the flip graph of regular triangulations and perform 
flips as in Section 8.3.1 of Ref. DLRS10 . This is the analogue of a Grobner walk 
in the setting of secondary fans. 

To find the secondary cone in the link of u given by a ray v, we compute the 
subdivision induced by u + sv for sufficiently small e > 0. Such a vector u + ev is 
represented symbolically in a way similar to a matrix term order in Grobner basis 
theory. 
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3. Resultants with specialized coefficients 

For some applications such as implicitization we need to compute resultant va- 
rieties while specializing some of the coefficients to constants. This problem was 
studied in Refs. |EKP07[ IEFKP11] for the case when the resultant variety is a 
hypersurface. In that case, the Newton polytope of the specialized resultant is 
the projection of the resultant polytope, and the authors computed the projection 
of resultant polytopes using Sturmfels' formula for vertices of resultant polytopes 
[Stu94[ Theorem 2.1] and beneath-bcyond or gift- wrapping methods for computing 
convex hulls. In our language, computing a projection of a polytope is equivalent 
to computing the restriction of the normal fan to a subspacc. 

In tropical geometry, specialization of certain coefficients amounts to taking sta- 
ble intersection of the tropical resultant with certain coordinate hyperplancs. In 
this section we first define the specialized tropical resultants and then present an 
algorithm for their computation. 

A polyhedral complex in R™ is called locally balanced if it is pure dimensional 
and the link of every codimension one face positively spans a linear subspace of R n . 

Definition 3.1. Let T\ and T 2 be locally balanced fans in R™. We define the 
stable intersection as the fan 

T\ n st T 2 := {Ci n C 2 :(C*i, C 2 ) 6Jix/ 2 and 

supp(link Cl (J- 1 )) - supp (link C2 (.F 2 )) = R™} 

with support 

supp(J- 1 n st T 2 ) = {uj e R" : supp(link w (J-i)) - supp(link w (J- 2 )) = R"}. 

If in addition T\ and T 2 are balanced then the stable intersection inherits multi- 
plicities from link tJ (J : i) and link^J-^) as follows: 

mult^-Fi (~l st T<i)~ 
mult Cl (link w T\ ) ■ multc 2 (link^ T 2 ) ■ [Z n : (Z n n RC*i ) + (Z" n RC 2 )] 

Ci,C 2 

where the sum runs over C\ G link (i ,(J r i) and C 2 £ link^J-^) such that a/ 6 C\ — C 2 
for a fixed generic vector uj' S R". 

Notice that the support of T\ n st J- 2 depends only on supp(J r i) and supp(J r 2 ). 
We will therefore extend the definition of stable intersections to intersections of 
supports of locally balanced fans and regard them as subsets of R" . 

For proofs of the following six statements, of which some are known to the 
community already, we refer to the upcoming paper Ref. [JYllj . 

Orthogonally projecting a polytope onto a linear space is equivalent to stably 
intersecting the tropical hypersurface of the polytope with the linear space. 

Theorem 3.2. Let P C R™ be a polytope, L C R™ be a linear subspace, and 
tv : R n — > L be the orthogonal projection. Then 

T(ir(P)) = (T{P) n st L) + L x . 

Lemma 3.3. For any locally balanced fans T\, T 2 , andJ-3, we have 

(1) n st T 2 ) n st T 3 = Tx n st (T 2 n st J3) 

(2) (supp(J r i) u supp(j" 2 )) n st suppp^) = supp^i n st T3) u supp(J" 2 n st -F3). 
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Lemma 3.4. For locally balanced fans T\ and T2 in R n 

supp(J" 1 n 5 * T 2 ) = (J Ci n C 2 . 

dim(Ci+C 2 )=n 

Corollary 3.5. Lef ^4 and B be locally balanced polyhedral fans in R™. Then 

link*, (A) n st Hnk w (S) = link w (A n flt B). 

Proposition 3.6. TTie stable intersection of two locally balanced fans is either 
empty or a locally balanced fan whose codimension is the sum of the codimensions. 

Lemma 3.7. Let I be an ideal rok[ii,i2,. . . ,x n ]. Then 

supp(T(/)) n st {x : X! = 0} = supp(T«/) + (xi - a))) 

where (I) is the ideal in \s.(a)[xi, X2, ■ ■ ■ ,x n ] generated by I. 

Definition 3.8. Let S = (Si, . . . , Sk) with Si C {1, . . . , m,} represent a choice of 
points in the configuration A. The coefficients of the monomials indexed by S are 
called specialized. Let Ui := {x <E l m ' : Vj 6 5* : Xj = 0} and U s ■= Ui=i u i- We 
define the specialized tropical resultant TTZs(A) := TTZ(A) P\ s t {Us}- 

We will use the following proposition to justify the word "specialized" : 

Proposition 3.9. Let A and S be as in Definition \3.8[ Let L be the ideal of 
1Z(A) and add to it, to obtain an ideal J, for each specialized coefficient Cj the 
binomial Cj — jj where jj is a parameter. We define the specialized resultant variety 

lZs{A) := V(J) C Jli=i(^*)" ii where K is the field of rational functions in the jj 
with coefficients in C. Then the tropicalization of lZs(A) is TH-s{A). 

Proof. The statement follows from Lemmas 13.71 and 13. 3f lV □ 

The computation of the tropicalization of 7^.5 (A) can be performed using Buch- 
bcrger's algorithm as explained in Ref. [BJS + 07] over the field of rational functions 
in the 7's. During this computation finitely many polynomials in the 7's appear 
as numerators and denominators of the coefficients. Substituting constant values 
for the 7's will give the same computation unless one of these polynomials van- 
ish. Hence specializing 7's to values outside a hypersurfacc in (C*) s will lead to a 
specialized tropical resultant variety. This explains the word "specialized" . 

If TTZs(A) is nonempty, then its codimension can be computed using Proposi- 
tion [313 and the codimension formulas from Section [2.41 Thus it remains to give an 
algorithm for checking if the specialized resultant is empty. Recall that m := J^. 
is the total number of points in A. 

Lemma 3.10. Let A and S be as in Definition \3.8[ Define the extended tuple 
B = . . . ,-Bfc) where Bi consists of bi^ = (di,j,Vi j) G Z n x Z m_ ' s ', with Vij £ 
j/n—\S\ b e i n g if j g Si and a standard basis vector otherwise. If the standard 
vector is chosen differently for every non- specialized coefficient then 

TH S (A) ^ m(B) = R rn . 

Proof. According to Lemma [3T4l TTZs(A) ^ if and only if there exists a cone C C 
TTZ(A) such that Us + C = K m where Us is as defined in Definition 13. 81 According 
to the simple description of tropical resultants in Theorem l2.9l we may assume that 
C has the form IR> {ey : a%j Ei} + rowspace(Cay(„4)). Equivalently, the stable 
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intersection is nonempty if and only if there exists a choice E such that M.>o{eij : 
fly ^ Ei} +rowspace(Cay(„4)) + Us has dimension m. Applying Theorem [279] to B, 
this is equivalent to T1Z(B) being full-dimensional, since rowspace(Cay(^4)) + Us = 
rowspace(Cay(£>)). □ 

Combining Lemma 13.101 and the results from Section 12.41 about codimcnsion 
computations, we get a polynomial time algorithm for deciding if a specialized 
result is nonempty. Another consequence of the lemma is the following algorithm 
for checking membership of a point in a specialized tropical resultant. 

Algorithm 3.11. (SpecializedResultantContains(A, S,lo)) 

Input: A tuple A of point configurations and a choice S of specialized coefficients. 
A vector wCP. 

Output: "True" if uj eTU s (A), "False" otherwise. 

• Compute the mixed subdivision of A induced by uj by computing the regular 
subdivision o/Cay(„4) induced byuj. 

• For each fully mixed cell: 

— construct a subconfiguration A! of points involved in the cell. 

— Return "True " if the specialized resultant of A' is nonempty. 

• Return "False". 

Proof. By Lcmma l3~37 2). Proposition ^ . 1 71 and Corollarv l3.5l we have that the sup- 
port of link w (7~7?.s(-^)) is the union of supports of T1Zs(A'), under the appropriate 
identification of TTZs{A') as a subset of M. m , where A' runs over all fully mixed 
cells of the mixed subdivision of A induced by u>. Hence u> 6 TlZs(A) if and only 
if one of TlZs{A') is nonempty. □ 

Algorithm 3.12. (NonTrivialVectorInSpecializedResultant(A, S ')) 

Input: A tuple A of configurations, a choice S of specialized coefficients such that 

U s n rowspace(Cay(.4)) C TK S (A). 

Output: A vector uj G TlZs{A) \ rowspace(Cay(^4)) 

• For each E = (Ei,En, . . . ,Ek) : Ei is a two- element subset of Ai, 

— Let C = K> {eij : i £ Ej} + rowspace(Cay(^l)). 

— J/dim(C + Us) = n and Us n C ^ Us H rowspacc(Cay(^4)) then 

* Find among the generators of Us H C a vector v outside the 
subspace Us H rowspace(Cay(^l)). 

* Return v. 

The following recursive algorithm finds a perturbed point in a starting cone for 
the specialized tropical resultant TlZs(A). 

Algorithm 3.13. (StartingPoint(A, S)) 

Input: A tuple A of configurations, a choice S of specialized coefficients such that 
TTls(A) + 0. 

Output: A vector w e £ Q(e) m such that for every fan structure ofl~lZs{A) defined 
over Q it holds that for e > sufficiently small, u> £ is in a facet of TlZs(A). 

• If dim(T!Zs(A)) = dim(Us H rowspacc(Cay(^l))) ; then return b\ + eb 2 + 
■ ■ ■ + e t_1 &t where b\, . ■ . , b t is some basis of Us n rowspace(Cay(^l)). 

• Compute an oj £ T1Zs{A) \ rowspace(Cay(^l)) using Alaorithm \3. 121 

• Compute the subdivision A w o/Cay(„4) induced byui. 

• For every fully mixed cell in A w . 
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— Let A' be the subconfiguration of the involved points. 

— Let S' be the restriction of S to A' . 

— Lf codimcnsion(T7?.s' (A')) = codimcnsion(T7?.s(.4,)) then 

* Return tu + e ■ StartingPoint(.4', S'). 

Proof. The correctness of the algorithm follows from the facts that the link of u> 
in the tropical resultant is the union of tropical resultants corresponding to the 
fully mixed cells in A w (Proposition I2.17|) , and that taking links commutes with 
taking stable intersections ( Corollary I3.5j) . and because the returned value from the 
recursive call is, after expansion with zeros, a vector outside of the secondary cone 
of w. □ 

We now turn to the problem of enumerating all facets in TTls(A) considered 
as a subfan of the restriction of the secondary fan of Cay(.4) to the subspace Us- 
While stable intersections are not in general connected in codimension 1, this is 
indeed the case f or TlZs(A) since it is defined by a prime ideal. See the argument 
m Ref. |BJS+07] using Kleinman's version of Bcrtini's Theorem. 

The output of Algorithm 13. 131 can be converted into a secondary cone in TTZ(A) 
containing uj £ in its relative interior, for example by computing a maximal secondary 
cone containing w e and taking the face containing uj e in its relative interior. 

The polynomial time algorithm of Theorem 12.291 for computing codimension a 
resultant varieties can be used for finding generic points in TTZ(A) in polynomial 
time. Simply remove points from A as long as possible without dropping the 
dimension. When done we have exactly two points left from each configuration of 
A. We then compute a generic point in TTZ(A) using Theorem l2.91 possibly using a 
symbolic e. It is unclear if a polynomial time algorithm exists for finding a generic 
point in specialized tropical resultants. 

Following the approach of Ref. jJenlOj discussed in Section 12.51 we are left with 
the problem of computing the link of a ridge in TTls{A). If the subspace Us 
had been generic enough to intersect the lineality space of T1Z(A) transversely, i.e. 
codim([/g nrowspace(Cay(^4)) = codim([/g) + codim(rowspace(Cay(.A))), then the 
link would be combinatorially equivalent to the link in T7l(A) and the support of 
the link would be a union of resultant fans of subconfigurations (Proposition 12 . 1 7j) 
where each fan can be found using Theorem 12.91 If Us is not generic, then com- 
puting a stable intersection with U s is required for finding the link in TlZs(A) 
(Corollary [33]). This is Algorithm 13. 141 below. Another approach is to compute the 
restriction of the secondary fan of each fully mixed subconfiguration to Us- We then 
get the resultant fan as certain rays of the secondary fan. This is Algorithm 13.151 

Algorithm 3.14. StableLinkfA, S,uj) 

Input: A tuple A of configurations, a choice S of specialized coefficients, a vector 
uj £ K" in the relative interior of a ridge R of TlZs(A). 
Output: A vector in each facet o/link w (7~7^s(^l)). 

• Let d be the dimension of TTZ(A) D s t {Us}- 

• Compute the subdivision A w of Cay(„4) induced by oj. 

• Z:=0. 

• For every fully mixed cell in A w 

— Let A' be the subconfiguration of involved points in the cell. 

— For each E = {E±, E2, ■ ■ ■ , E^) : E± is a two-element subset of A[, 

* Let C = R>o{ei.j : i £ Ej} + rowspace(Cay(^4)). 
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* If dim(Us + C) = m and dim(i7s fl C) = d then 

■ Let V be a set of one or two vectors in Us H C such that 
(Us f~l C) + span(i?) is positively spanned by V U span(-R). 

■ l:=lUV 

• Return I. 

Algorithm 3.15. StableLinkfA, S,u) 

Input: A tuple A of configurations, a choice S of specialized coefficients, a vector 
to £ M. n in the relative interior of a ridge R ofl~TZs(A). 
Output: A vector in each facet of]mk U! (T'R's(A)). 

• Let d be the dimension ofTTZ(A) r\ st {Us}- 

• Compute the subdivision A u of Cay (A) induced by u>. 

• I :=0. 

• For every fully mixed cell in A w 

— Let A' be the subconfiguration of the involved points of the cell. 

— If the codimension of the lineality space of the restriction IF of the 
secondary fan £(Cay(.4')) to Us is m — d, then 

* Choose v such that v completes span(i?)n£/s to a generating set 
of the lineality space of T . 

* If SpecializedResultantContains (A' , S,r) then I := I U {v, —v}. 

— else 

* Compute all maximal cones in J- (by traversal) . 

* For each ray r in J- , if SpecializedResultantContains (A' , S,r) 
then 1:=IU {r}. 

• Return I. 

The above algorithm is to be read with proper identifications - when restricting 
to A' the vectors in K m need to be truncated accordingly, and so does the set S, 
and r needs to be expanded when adding it to I. When adding vectors to I, it is 
advantageous to choose the vectors as primitive vectors orthogonal to the span of 
the ridge so that duplicates can be removed easily. 

If Us is high dimensional, a typical situation is that each subconfiguration is a 
number of edges and a triangle. In this case there are only few choices E to run 
through in Algorithm 13.141 For lower dimensional Us there can be many choices 
of E but with many of the contributions to the stable intersection being the same. 
See Example 13.161 In such a case Algorithm 13.151 performs better. In general 
it is difficult to predict which algorithms is best. In our implementation we use 
Algorithm 13.151 and Algorithm 13.141 only when there is no specialization. 

Example 3.16. Let A = (A 1 ,A 2 ,A 3 ) with 

A 1 = {(0,0), (0,1), (0,3), (1,0), (3,0)} 

A 2 = {(0,0), (0,1), (0,3), (1,0), (3,0)} 

A 3 = {(0, 0), (0, 1), (0, 2), (1, 0), (1, 3), (2, 0), (3, 1), (3, 3)}. 

Choosing the specialization S of every coefficient except the coefficient of the point 
(0,0) in each configuration, we get that TTZs(A) is a two dimensional fan with 
f-vector (1, 13, 17) living inside M. 3 C M 18 . The link at e n G M 18 consists of 4 rays. 
The traversal of TTZs (-A) takes 79 seconds if Algorithm 13.141 is used but only 5 
seconds if Algorithm 13.151 is used for computing the links. Algorithm 13.141 needs 
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to iterate through 2100 vertex pair choices at en, but much fewer for many of the 
other links. 

3.1. Implicitization using specialized resultants. In this section we will show 
that the tropicalization of a variety parameterized by polynomials with generic co- 
efficients can be computed using specialized tropical resultants. Let f\, /2, . . . , fk S 
C[a;f 1 , x^ 1 , . . . , x^ 1 ) be polynomials parameterizing a variety V in C k . Let T be 
the graph of the parameterizing map, defined by (yi — fi,y2 — fs, . ■ . ,yk — fk) in 
C[x ± 1 , xf \ . . . , x± 1 ,yi,y 2 , ■ ■ ■ , yk]- When /i, / 2 , ...,/* have generic coefficients, 
the tropical variety of V is the stable intersection of the tropical hypersurfaces of 
the polynomials y\ — /i , y% — fi , ■ • ■ , yk — fk- Since V is the closure of the pro- 
jection of r C C" x C k onto C k , by tropical elimination theory, we can compute 
the tropical variety of V as a projection of T(T). This approach was used in 
Rcfs. [STY071ISY08] . 

Another way to compute the tropical variety of V is by using specialized resul- 
tants. Let A ~ (Ai, A2, . . . , Ak) where Ai = supp(/;) U {0} for each i = 1, 2, . . . , k. 
Let S = (supp(/i), supp(/2), . . . , supp(/fc)) be a choice of points to specialize, and 
let Us be the subspace of Y\i=i x ^™ defined by setting the coordinates in S 
to 0. 

Proposition 3.17. With the notation above, T(V) = TlZs(A), i.e. the tropical- 
ization of a variety parameterized by polynomials with generic coefficients coincides 
with a specialized resultant. 

Proof. Let W be the incidence variety in rii=i( C *) A< x ( C *)" 

as seen m (1]), defined 

by equations of the form yi — gi where gi is a polynomial with the same support as 
fi but with indeterminate coefficients. Then the graph T is obtained by specializing 
the coefficients of gi to those of fi. Since the coefficients of fi were assumed to be 
generic, we get T{T) = T(W) n st U s . By tropical elimination, T(V) = T{T)+R n , 
and by the following lemma, this coincides with (T(W) + M") n s t Us, which is the 
specialized tropical resultant. Here we are using the fact that a linear projection 
of a fan is combinatorially equivalent to the Minkowski sum with the kernel of the 
projection. □ 

Lemma 3.18. Let J- be a locally balanced fan in W l . Let L and L' be linear 
subspaces o/K n such that L' C L. Then 

(Fn st L) + L' = (F + L r ) n st L 

In other words, stable intersection with a linear space commutes with Minkowski 
sum with a smaller linear space. 

Proof. Both (J 7 n s t L) + L' and (F + L') fl s t L arc empty if T + L has dimension less 
than n. Suppose this is not the case. Then both sets contain L' in their lincality 
space and consist of points of the form u + v <G K™ where u £ L' and v € J- n L are 
such that dim(link„(J r ) + L) = n. □ 

Since the tropical variety of the graph V only depends on the extreme monomials 
of the parameterizing polynomials, the following result follows immediately. 

Corollary 3.19. When using specialized resultants for implicitization, the extreme 
monomials of the input polynomials determine the tropical variety of the parame- 
terized variety, so we can safely disregard the non-extreme terms. 
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An advantage of using specialized resultant is that our description of special- 
ized tropical resultant is easier to compute than the stable intersection directly. 
Moreover, experiments show that the resultant description may speed up the re- 
construction of the Newton polytopc in some cases. See Section [5] for examples. 

Moreover, when the variety V is not a hypcrsurfacc, our resultant description 
gives a fan structure of T{V) derived from the restriction of a secondary fan to a 
linear subspace, which is the normal fan of a fiber polytopc. Tropical elimination 
does not give a fan structure for varieties of codimension more than one. 

3.2. Tropical elimination for specialized tropical resultants. As before, let 
A be a tuple point configurations in Z™ and S be the tuple of subsets to be special- 
ized. Let W be the incidence variety and TW be is tropicalization as in Section [2~T1 
Let Ws be a variety cut out by polynomials fi where the the coefficients of mono- 
mials in S have been specialized. Then fi, f%, . . • , fk may no longer form a tropical 
basis, but the tropicalization of Ws can be computed as the stable intersection of 
tropical hypersurfaces of /i, /2, ■ • ■ , /fc because the coefficients are assumed to be 
generic (or indeterminate). 

The specialized resultant is the projection of Ws onto the non-specialized coef- 
ficient variables, and we can compute this using tropical elimination theory, which 
gives the tropical variety as a union of cones. When the specialized tropical resul- 
tant is a tropical hypcrsurfacc, then we can reconstruct the normal fan of the dual 
Newton polytope using the methods in the next section. 

The tropical hypersurface of fi only depends on the Newton polytopc Pi of fi. 
The non-specialized points in Ai always contribute as vertices of Pi, but some 
specialized points of Ai may not. From this observation, we obtain the following 
result, which is not obvious from the resultant point of view. 

Lemma 3.20. If G Ai is a specialized point lying in the convex hull of other 
specialized points in Ai, then removing aij from Ai does not change the specialized 
tropical resultant because the Newton polytope and the tropical hypersurface of fi 
remain the same. 

In other words, we may disregard the non-vertices among the specialized points. 
Using this lemma, we may be able to reduce the amount of work for computing the 
specialized tropical resultant TlZs{A). 

4. Polytope reconstruction 

In this section we describe an algorithm for finding a fan structure on a tropical 
hypcrsurfacc T C R™ . Recall that the tropical hypersurface of a polytopc P C K™ 
is the set of ui € K™ for which there exist distinct p,q e P such that for any r € P, 
uj-p = uj-q<uj-r. In other words, the tropical hypersurface of a polytope is 
the union of the normal cones to the polytope at the edges. The multiplicity of 
a point in the relative interior of such a normal cone is the (lattice) length of the 
edge. The tropical hypersurface of a polynomial is the tropical hypersurface of its 
Newton polytope. 

The tropical hypersurface T will be presented to us as a finite collection of 
codimension 1 cones which may overlap badly but whose union will be T . What 
we wish to compute is a collection of codimension 1 cones such that the collection 
of all their faces is a polyhedral fan with support T. This fan is not unique unless 
we require it to be the coarsest - that is, that it is the normal fan of the polytope 
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defining T with its maximal cones removed. If the codimension 1 cones come with 
a multiplicity then an advantage of having the fan structure is that it is relatively 
easy to reconstruct the 1-skclcton of the polytopc defining T up to translation. 
Therefore we will consider the computations of a polytopc, its normal fan, and its 
tropical hypersurface to be equivalent in the following. 

One way to perform the polytope reconstruction is to use the beneath-beyond 
method for computing convex hulls. The key observation is that for any generic 
oj G R n the vertex face u (New(/)) can be computed using "ray shooting". See the 
papers Refs. [DFS07] and |CTY10j . The method we present in this paper uses the 
adjacency decomposition approach (see Section \2.5\i and the following algorithm 
for computing normal cones, at vertices, of the polytope defining T. 

Algorithm 4.1 (Rcgion(S,w)). 

Input: A collection S oj codimension 1 cones in R™ such that T := UcesC *s the 
support of a tropical hypersurface. A vector uj G R™ \ T. 
Output: The (open) connected component o/R\T containing to. 

• R :=R™. 

• For each C <E S : 

- While RHC ^ 0: 

* Find p G R n C. 

* Let h be the generic open half line from ui through p + X)"=i £le * ■ 

* Let A C S be the collection of cones with intersection point with 
h closest to uj, for e > sufficiently small. 

* Let for some D G A the half space H C R" be the connected 
component o/R" \ span(D) containing uj. 

* R:=RC\H. 

• Return R. 

Proof. The set R stays open and convex throughout the computation. At the 
end R PI T = 0. Each added constraint H for R is necessarily satisfied by the 
connected component because of its convexity. The symbolic perturbation of p and 
the convexity of R ensures that H is independent of the choice of D in A, as all 
cones in A must be parallel. In fact the set of constraints give an irrcdundant 
inequality description of the returned cone. □ 

Proposition 4.2. Let a be the number of facets of the closure of the returned cone 
of Algorithm \4-l\ The number of checks "RD C ^ " performed in algorithm is 
\S\ + a while the number of interior point computations "p G R D C " is a. 

Proof. The check is done for every cone in C G S. In addition, whenever the 
algorithm enters the body of the while loop, a facet constraint H is added to R, 
and an additional check U R H C ^ 0" and a computation of p is performed. □ 

The condition that the generic h intersects a given polyhedral cone C can be 
phrased as a condition on the ordering in which h intersects the defining hyperplancs 
of C . We can imagine moving a point starting from ui and along the half-line 
h, keeping track of which equations and inequalities defining C are satisfied and 
updating when a defining hypcrplane of C is crossed. Hence the implementation 
reduces to a check of the order in which h intersects two given hyperplancs. The 
perturbation in such a check is not difficult to handle symbolically. The check can 
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be used again to actually find a D in the algorithm with intersection point closest 
to 10. 

To apply the adjacency decomposition approach we must be able to compute 
a starting cone and move across codimcnsion 1 faces to find neighboring cones, 
while computing links of ridges is trivial for complete fans. To find a starting 
cone we guess a vector outside T and apply Algorithm 14.11 Suppose now that 
C is a full dimensional cone in the normal fan and u is a relative interior point 
on a facet of C with outer normal vector v. For e > sufficiently small, calling 
Algorithm 14. II with argument u + ev will give us the desired neighbouring cone. In 
our implementation we again use comparison of intersection points on line segments 
to find an e sufficiently small to avoid all hyperplanes appearing in the description 
of T. 

If we precompute generators for the cones in S then most of the checks for 
empty intersection with R can done without using linear programming, but rather 
for each defining hypcrplanc of R checking if the cone generators arc completely 
contained on the wrong side. In our current implementation the time spent on 
finding first intersection along the half-lines is comparable to the time spent on 
linear programming. We present two examples to illustrate the usability of the 
algorithm. 

Example 4.3. In Ref. |HSYY08] we computed the f-vector of the tropical hypcr- 
surface of the 2x2x2x2 hyperdetcrminant. Its support is the sum of a tropical 
linear space and a classical linear space in M 16 and is to easy write as a union of 
cones. We reconstruct the 25448 normal cones of the Newton polytope of the defin- 
ing equation in 163 minutes. Exploiting the 384 order symmetry as explained in 
Ref. [JenlOj we reduce the running time to 7 minutes for computing the 111 orbits 
of maximal cones. With suitable input files the following Gfan command |Jen| will 
compute the f-vector. Also see Section [5] 

andersSgureko : ~$ gf an_tropicalhypersurf acereconstruction -i troplinspc . f an 
— sum — symmetry <claslinspc_and_symmetry . txt I grep -Al F_VECT0R 
F_ VECTOR 

1 268 5012 39680 176604 495936 927244 1176976 1005946 555280 178780 25448 

Example 4.4. The implicitization challenge solved in Ref. jCTYlOj is to recon- 
struct the Newton polytope of the defining equation of a tropical variety given as a 
union of 6865824 cones. This 11-dimensional polytope lives in R 16 and has a sym- 
metry group of order 384. In Ref. [CTYlOj , a modified version of the ray-shooting 
method was used to produce coordinates of vertices at a rate of a few (2-5) min- 
utes per vertex. Each round took about 45 minutes found 10-20 vertices typically. 
However a lot more computation, with some human interaction and parallclization, 
over a period of a few months was required to make sure that all the vertices were 
discovered, and this was done by computing the tangent cone at each found ver- 
tex, up to symmetry. During the process most vertices were re-discovered multiple 
times. 

On this example our new implementation in Gfan spends approximately 1 minute 
for each call of Algorithm UTT} We estimate that the enumeration of the 44938 orbits 
would finish after 30 days of computation. With the new method, we do not need to 
process a vertex more than once, and we obtain all the facet directions as the rays 
in the normal fan and all the tangent cones as duals of the normal cones. Moreover, 
there is no post-processing needed to certify that all vertices have been found. 



26 



ANDERS JENSEN AND JOSEPHINE YU 



The method we just described does not make use of multiplicities. In fact, it is 
not necessary that the fan is polytopal, or even locally balanced. We only require 
that each connected component of the complement of T is convex. 

Before settling with Algorithm 14. II we also experimented with storing the codi- 
mension one cones in a binary space partitioning tree (BSP tree). The tree would 
be built at initialization, and the connected components of the complement could 
be computed by gathering convex regions stored in the tree. This method worked 
as well as Algorithm 14.11 in small dimensions, but in higher dimensions, like the 
examples above, Algorithm 14.11 would always perform better. In Example 14.31 the 
difference would be a factor five without exploiting symmetry. But in Example 14.41 
the number of required nodes of the tree would grow too large to have any chance 
of fitting in memory. The intuition behind the explosion in complexity is that cones 
(for example, simplicial cones of codimension one) in a higher dimensional space 
have larger chances of intersecting a fixed hypcrplane. Therefore in the process of 
building the BSP tree, a codimension one cone from the input will meet many other 
hyperplanes coming from other cones, causing an explosion in the number of nodes 
in the BSP tree. 

5. Comparison of algorithms 

In this section, we consider various algorithms and compare the combinatorial 
complexity of the output (e.g. f-vector) and running time (recorded on a laptop 
computer with a 2.66 GHz Intel Core i5 processor and 8GB of memory). All 
implementations are single threaded, done in C++ using eddlib Fuk05j and SoPlex 
|Wun96j . and will be part of Gfan in its next release, unless otherwise noted. The 
combinatorial complexity of the output is essential for a fair comparison since not 
the same amount of effort went into making each of the implementations fast. A 
lot of effort went into the implementation of Algorithm 14.11 and the secondary fan 
computation because of their broad range of applications, while less optimization 
effort has gone into algorithms specific to tropical resultants. 

In general, the software Gfan uses the max convention for tropical varieties and 
Grobner fans. However, for the fact that the secondary fan of a point configuration 
is a coarsening of the Grobner fan of the associated binomial (lattice) ideal, we need 
the subdivisions to be defined with respect to min if the initial ideals are defined 
with respect to max. Therefore Gfan uses min for secondary fans. As tropical 
resultants are subfans of secondary fans, we chose to use min in this paper for 
tropical addition. 

Hypersurfaces. Let us first consider the case where the resultant variety 1Z(A) is 
a hypcrsurface. Following is a list of different methods for computing the resultant 
polytope (or its tropical hypersurface or its normal fan). 

(1) Enumerating the vertices of the secondary polytope of Cay(^4), and then 
using Sturmfels' formula |Stu94[ Theorem 2.1] to obtain the vertices of 
the resultant polytope. For our timings we only do the first part of the 
computation using the Gfan command 

gf an_secondaryf an <cayley.txt 

(2) Computing the tropical hypersurface of the resultant as a subfan of the 
secondary fan by fan traversal using the methods described in Section [2.51 

gf an_resultantf an — vectorinput <tuple.txt 
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(3) Constructing the normal fan of the resultant polytope from the simple 
description of the tropical resultant as a union of cones as in Theorem l2.9l 
Our implementation in Gfan uses Algorithm 14.11 for this. 

gf an_resultantf an — vectorinput — projection <tuple.txt 

(4) For a generic direction using Sturmfels' formula |Stu94[ Theorem 2.1] for 
finding the optimal vertex of the resultant polytope in that direction and 
combining this approach with the beneath-beyond convex hull algorithm 
for recovering the whole polytope. This method has recently been imple- 
mented. See Ref. |EFKPli] . Unfortunately the current interface for the 
implementation only handles implicitization type problems and cannot han- 
dle Examples (a) through (f) below. We intend to make comparisons with 
it for these examples in the near future. 

For the third method, one can also use other methods for reconstructing a poly- 
tope from its tropical hypersurfacc, such as ray-shooting/beneath-beyond and BSP 
trees, as discussed in Section^ although we found Algorithm l4.1l to perform better, 
especially for polytopcs of dimension 5 or more (compared to beneath-beyond in 
iB4e |Hug06| and BSP). 

For Example 12.101 above, each of the first three methods finished in under one 
second in Gfan. We will present more challenging examples below. 

Example (a). 

/ /0 1 3\ /0 l\ /0 2 2\ /l 2 2 

.4= 1,0 2 1,2 1 2 , 2 3 

\\1 1 1/ \3 2 0/ \3 1 1/ \1 2 



Method/fan 


F-vector of output 


Timing 


(1) secondary fan 

(2) traversing tropical resultant 

(3) normal fan from simple description 


1 10432 55277 106216 88509 27140 
1 5152 21406 28777 12614 
1 78 348 570 391 93 


467 s 
733 s 
1.4 s 



Example (b). 

. _ ( (0 1 3Wl 2 3 3\/0 1 2 3\\ 

v v° 1 2 ^'v 1 2 v'v 1 1 3 // 



Method/fan 


F-vector of output 


Timing 


(1) secondary fan 

(2) tropical resultant 

(3) normal fan 


1 3048 38348 178426 407991 494017 304433 75283 
1 2324 26316 106083 197576 173689 58451 
1 56 497 1779 3191 3018 1412 249 


506 s 
1238 s 
6 s 



Example (c). 
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2 




2 


3 


3 




3 


2 


2 






3 1 


2 


3 





1 


3 




2 


3 










1 







1 


1 


2 




V 


\l 2 


V 






3 








3 


<V 






3 


-V 









V 


) 



Method/fan 


F-vector of output 


Timing 


(3) normal fan from simple descr. 


1 937 5257 11288 11572 5589 985 


55 s 



In Example (c) we were not able to compute the secondary fan and the resultant fan 
with the secondary fan structure due to integer overflow in intermediate polyhedral 
computations. Gfan has been designed to work well for Grobner fans, where the 
degrees of the polynomials is never very large, since that would prevent us from 
computing a single Grobner basis anyway (except for binomial ideals) . In Example 
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(c), a primitive normal vector of a codimension 1 cone of the normal fan of the re- 
sultant is (-32, 0, 32, 27, 0, -27, 25, -25, 0, 0, 51, -51, -87, 0, 87), showing that the 
resultant has degree at least 32+27+25+51+87. On such examples overflows typi- 
cally arise when trying to convert an exactly computed rational generator of a ray 
to a primitive vector of 32 bit integers. Algorithm 14.11 will show similar behavior 
on other examples, for example when converting u p" to a vector of 32 bit integers. 
We intend to fix these implementation problems in the future. 

Hypersurfaces with Specialization. If the specialized resultant is a hypersur- 
face, then we can compute its tropical variety using the following methods. 

(1) Compute TlZs(A) as a subfan of the restriction secondary fan to a subspace 
Us by fan traversal using the algorithms in Section |3] 

gf an_resultaiitf an — vectorinput — special <tuple_and_spcvec . txt 

(2) Compute the stable intersection TlZs(A) = T1Z(A) fl s t {Us} as a union of 
cones, using the simple description from Theorem l2.9l and the characteriza- 
tion of stable intersections from Lemma 13.41 Then reconstruct the normal 
fan of the dual polytope using Algorithm 14. II 

gf an_resultantf an — vectorinput — special — projection <tup_and_sv.txt 

(3) Compute the specialized tropical resultant as a union of cones using sta- 
ble intersection of hypersurfaces and tropical elimination theory as in Sec- 
tion 13.21 and reconstruct the normal fan of the dual polytope using Algo- 
rithm dJ] We combine the commands (see also Ref. |SY08] V. 

gf an_tropicalstartingcone — stable >startingcone .txt 

gf an_tropicaltraverse — stable <startingcone . txt >stable.fan 

gf an_tropicalhypersurf acereconstruction — sum -i stable. fan <lnspc.txt 

(4) For a generic direction, Sturmfels' formula |Stu94[ Theorem 2.1] gives the 
optimal vertex of the resultant polytope in that direction, which can then be 
projected to get a point in the Newton polytope of the specialized resultant 
polynomial. This can be combined with the beneath-beyond convex hull 
algorithm for recovering the whole polytope. This method has recently been 
implemented as described in Ref. [EFKPli] , This is the implementation 
used in the timings below. 

In Ref. jEKP07j . the authors proposed computing a silhouette or a projection of 
the secondary polytope. This is dual to computing the restriction of the secondary 
fan to a subspace. We provide the results and timings of this dual computation for 
comparison. 

In the following examples the non-black columns are specialized. 
Example (d). 

A-ff° 1 ^ f° 1 1 2 \ f° 1 1 2 ^ 
V \° 1 V ' V 1 V ' \° 1 2 V / 



Method/fan 


F-vector 


Timing 


Restriction of secondary fan 

(1) traversing tropical resultant 

(2) normal fan from stable intersection 

(3) normal fan from tropical elimination 


1 372 2514 5829 5661 1976 
1 126 476 561 212 
1 25 127 250 211 65 
1 25 127 250 211 65 


26 s 
14 s 
0.7 s 
1.4 s 
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Example (e). 

A-((° 1 !W° 1 1 2W0 1 1 2\\ 
V \° 1 V ' V 1 1 V ' V° 1 2 V / 



Method/fan 


F-vector 


Timing 


Restriction of secondary fan 

(1) traversing tropical resultant 

(2) normal fan from stbl. inters. 

(3) normal fan from trop. elim. 


1 709 6955 24354 39464 30226 8870 
1 469 3993 11296 12853 5040 
1 29 209 597 792 485 110 
1 29 209 597 792 485 110 


116 s 
320 s 
1.3 s 
3.2 s 



Example (f). 



//l 1 2 3\ /() 1 1\ (1 1 2 3\ (I 1 3 3\\ 
A = 2 2 3 2,1 2 1 1,1 3 3 2,0 2 1 

\\0 2 1 2) \0 2 1 3/ \1 1 1/ \3 2 1 1/ J 



Method 


F-vector 


Timing 


(2) 
(3) 


1 1566 19510 98143 265202 424620 413455 238425 73741 9156 
1 1566 19510 98143 265202 424620 413455 238425 73741 9156 


798 s 
974 s 



Implicitization of hypersurfaces. Implicitization is a special case of the spe- 
cialized resultants, and we compare the three methods as before. 
Example (g). (Implicitization of a bicubic surface. Example 3.4 of Ref. |EK05j ) 

. _ f [0 1 2 3W0 1 2 3\ 
\\0 1 2 3 o)'\0 1 3 Q)' 



(0 001111222233 3\\ 
^0 1201230123123// 



Method/fan 


F-vector 


Timing 


No interior points 


Restriction of secondary fan 


1 26 66 42 


5 s 


2 s 


(1) traversing tropical resultant 


1 13 17 


16 s 


4 s 


(2) normal fan from stable inters. 


15 9 6 


171 s 


9 s 


(3) normal fan from tropical elim. 


15 9 6 


0.4 s 


0.4 s 


(4) beneath-beyond 


15 9 6 


< 0.1 s 


< 0.1 s 



As we saw in Corollarv l3.19l removing the non-extreme monomials from the param- 
eterizing polynomials does not change the resultant polytope, and in this example, 
this also does not change the restriction of the secondary fan. However, doing so 
speeds up the computations, as seen on the right most column. 
Example (h). (Implicitization of a hypersurface in four dimensions) 

//0 2 4\ /0 1 2 3\ /0 2 3 4\ /0 4 4\\ 
A= 2 4 1,0 2 2 0,0 4 1,0 2 2 3 

\ \0 2 4 1/ \0 1 4 lj \0 2 4 2/ \0 4 2 3/ / 



Method/fan 


F-vector 


Timing 


(2) normal fan from stable intersection 

(3) normal fan from tropical elimination 

(4) beneath-beyond 


1 111 358 368 121 
1 111 358 368 121 
1 111 358 368 121 


9 s 
2.6 s 
1.3 s 



For (3), computing the polytope from the tropical hypersurface using ray-shooting 
and beneath-beyond took 47 s in the Trim implementation |SY08| using the library 
iB4e |Hug06| on a slightly slower machine. 

Example (i). (Implicitization of a hypersurface in five dimensions) 
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.4: 





f° 


1 


3 







1 


4 







2 


2 


I 


\fl 


2 


■1 



/0 1 3\ /0 2 3\ /0 1 2 3\ /0 2 4\\ 
2 3 1 4 2 1 4 2 4 1 3 
0113'0111'0010'0343 

Vo 1 2 3/ \0 4 2 3/ \° 1 3 3 / \° 1 3 W 



Method/fan 


F-vector 


Timing 


(2) normal fan from stable inters. 

(3) normal fan from tropical elim. 

(4) beneath-beyond 


1 5932 23850 35116 22289 5093 
1 5932 23850 35116 22289 5093 
1 5932 23850 35116 22289 5093 


351 s 
184 s 
1241 s 



For (3), timing includes 17 seconds for computing the specialized tropical incidence 
variety. The normal fan reconstruction computation in Trim with iB4e took 3375 
seconds on a slightly slower machine. 

Non-hypersurfaces. When 1Z(A) is not a hypcrsurface, the only method we know 
for computing T1Z(A) with a fan structure without knowing the defining ideal is 
to traverse the secondary fan of Cay(„4) and enumerating just the secondary cones 
whose RMS contains a fully mixed cell. There are other descriptions of tropical 
resultants as a set, such as Theorem 12. 91 but none gives a fan structure. 
Example (j). 



A = 



2 4 
4 1 1 



3 5 5 
1 4 



3 4 5 
1 5 2 



1 2 

4 3 5 



Method/fan 


F-vector 


Timing 


Secondary fan 
Traversing tropical result. 


1 8876 72744 222108 322303 225040 60977 
1 968 4495 6523 3000 


478 s 
81 s 



We used, respectively, the commands: 

gf an_secondaryf an <cayley.txt 

gf an_resultantf an — vectorinput <tuple.txt 

Non-hypersurfaces with Specialization. The only method here is to traverse 
T1Zs{A) as a subfan of a restriction of the secondary fan using the algorithms in 
Section [3J 
Example (k). 



A 



2 4 
4 1 1 



3 5 5 
1 4 



3 4 
1 



5 

5 2 



12 

4 3 5 



Method/fan 


F-vector 


Timing 


Restriction of secondary fan 
Traversing spec, tropical result. 


1 4257 23969 48507 42260 13467 
1 310 831 533 


256 s 
81 s 



We used, respectively, the commands: 



gf an_secondaryf an 
gf an_resultantf an 



-restrictingf an subspace.fan <cayley.txt 
-vectorinput — special <tup_and_sv.txt 



5.1. Conclusion. The new method of using adjacency decomposition with Algo- 
rithm l4TT1 for constructing the normal fan of a polytopc from it tropical hypcrsurface 
works very well in practice. Our implementation of it is much faster than any exist- 
ing implementation of the beneath-beyond method with ray-shooting, and we think 
the gap will widen even more in higher dimension since this new method scales well 
multi-linearly with respect to the number of cones in input and the number of 
vertices and edges of the output polytope, as shown in Proposition ^. 21 
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The normal fan reconstruction method can be used together with either the 
simple description of tropical resultants (Theorem l2.9[) or tropical elimination (Sec- 
tion [22]) f° r computing resultant polytopes efficiently. Traversing the (specialized) 
tropical resultant as a subfan of (a restriction of) the secondary fan of the Cayley 
configuration is combinatorially interesting but not computationally competitive. 

For implicitization, the beneath-beyond method from Rcf. jEFKPli] works faster 
than any of our "tropical" methods when output polytope is low dimensional, while 
our methods seem to have an advantage in higher dimension (5 or more). However, 
the method of Ref. [EFKPlT] may have an advantage when there are many spe- 
cialized points in the input configurations, as the number of cones in the tropical 
description increases rapidly. See the last problem in Section [6] below. 

For resultant varieties of codimension higher than one, whether specialized or 
not, we only know of one method for computing the tropicalization as a fan, without 
knowing the defining polynomials, which is to traverse the secondary fan of the 
Cayley configuration or a restriction of it to a subspace. 

6. Open problems 

Combinatorial classification of resultant polytopes: For 1-dimensional 
point configurations, the combinatorics of the resultant polytope only de- 
pend on the (partial) order of the (not necessarily distinct) points in each Ai 
[GKZ94] , so a combinatorial classification is easy to obtain. No such classi- 
fication is known even for point configurations in Z 2 . A concrete problem is 
to classify 4 dimensional resultant polytopes combinatorially. This 
was done for 3 dimensional resultant polytopes by Sturmfels |Stu94| . and 
only one dimensional point configurations were needed for this case. To 
understand the 4 dimensional resultant polytopes, we need to work with 
the case A = (A\, A2, A3) where each Ai consists of three points in Z 2 that 
are not necessarily distinct. How can we stratify the space of tuples A's 
according to the combinatorial type of the resultant polytope? 

Finding a point in the specialized tropical resultant: Is there a poly- 
nomial time algorithm for finding a generic vector lo € Q(e) m in the special- 
ized tropical resultant? For non-specialized tropical resultants, the polyno- 
mial time algorithm for codimension from Section 12.41 can also be used to 
find a generic point, by Theorem 12.91 

Improved description of specialized tropical resultants: Combining the 
descriptions of tropical resultants in Theorem 12.91 and stable intersections 
in Lemma 13^41 we get a specialized tropical resultant as a union of cones. In 
computations, we need to go through a list of JliLi ("2) choices of tuples of 
pairs from A j , many of which do not contribute to a facet of specialized trop- 
ical resultant. Give a combinatorial characterization for the choices 
of the tuples of pairs that contribute to a facet. Corollary 13. 191 and 
Lemma 13.201 are results in this direction. 
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